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Abstract 

In this paper excursions of a stationary diffusion in stationary state 
are studied. In particular, we compute the joint distribution of the 
occupation times I t and J t ^ above and below, respectively, the 
observed level at time t during an excursion. We consider also the 
starting time gt and the ending time dt of the excursion (straddling 
t) and discuss their relations to the Levy measure of the inverse local 
time. It is seen that the pairs ') and (t — gt,dt — t) are 

identically distributed. Moreover, conditionally on + l\ ' = v, 
the variables I t and l\ ^ are uniformly distributed on (0,v). Using 
the theory of the Palm measures, we derive an analoguous result for 
excursion bridges. 

MSC: 60J60; 60G10 
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1 Introduction 

In the literature there are many examples of cases such that the occupation 
time for one diffusion is identical in law with the first hitting time for another 
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diffusion. For these see, e.g., [Zl El H21 EDI EH El EH ESI EIJ - The last one of 
these references contains a survey of known identities. 

In this paper we prove a new kind of identity between hitting and occu- 
pation times for stationary diffusions in stationary state. An example of such 
a diffusion is reflecting Brownian motion on K+ with drift — /i < 0. As is well 
known, this process is stationary having the exponential distribution with 
parameter 2/j, as its stationary probability distribution. We take the whole 
of R to be the time axis and use, for a moment, the notation {X t : t G IR} 
for this process. For fixed t G IR let 

9t := sup{s < t : X s = 0}, d t : = inf{s > t : X s = 0}, (1.1) 

and 

rdt pdt 

1+ :— / lsx^xAds, := / li n ^x.<xA ds. 



h '■= / J -{x s >x t }«6, J-t ■= / ±{o<x s <x t } 
' gt Jgt 

Then from [29j it can be deduced that 

(li + \lt ] ) = (t-g t ,d t -t). (1.2) 

For an explicit form of this distribution and its Laplace transform, see Section 
5, Example 5.1. The main result of the present paper states that ()1.2|) is valid 
for all stationary diffusions in stationary state. By time reversal, the variables 
, t — gt, and d t — t are identically distributed. Our proof of (|1.2|) is 
purely computational and does not, unfortunately, provide any probabilistic 
explanation of the identity. In the case of the reflecting Brownian motion 
with drift we have an alternative, more probabilistic, proof based on the Ray- 
Knight theorems. It is fairly easy to deduce also in the general case using 
the symmetry and time reversal as is done in [2H1, that the expectations of 
If , t — gt, and d t — t are equal. 
It is of interest to recall the identity due to [3] and [12] because the 
variable d t — t can also be found in this identity. Let {Bf^ : t > 0} be a 
Brownian motion with drift fi > starting from and 

H X {B M ) := M{t : B^ = x} 

the first hitting time of the level x > 0. Then the Biane-Imhof identity states 
that 

l {B (, )<0} ds ( = } H X (B^), (1.3) 
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where A is exponentially (with parameter 2/i) distributed random variable 
independent of . By spatial homogeneity, we have 



since the distribution of X t is exponential with parameter 2/i. 

From the identity (|1.2jl using the theory of Palm measures applied for 
excursions we derive a new interesting result for excursion bridges. To explain 
this, let X(°'''°) be an excursion bridge of length I for excursions from to 
for an arbitrary diffusion X. Further, let U be a random variable having the 
uniform distribution on (0,/) and assume that U is independent of X^°' l '°\ 
Then the occupation times 



are identically uniformly distributed on (0,/). For a standard Brownian ex- 
cursion, i.e., for the 3-dimensional Bessel bridge, this fact can also be ex- 
plained via Vervaat's path transformation. 

The paper is organized so that there are two main sections and two shorter 
section with remarks and examples. The first main section is devoted to the 
identity ()1.2|) and is divided into four subsections. In the first subsection some 
preliminaries on stationary diffusions are given. After this the joint distri- 
bution of go and do is derived. In particular, it is seen that this distribution 
belongs to a special class of two-dimensional distributions characterized by 
the fact that the density (when it exists) is a function of the sum of the 
arguments only. The proof of the main identity (|1.2|) is presented in Section 
12.31 After this, in Section EH we give an alternative proof of (II. 2j) for RBM 
with drift based on the Ray-Knight theorems. In the second main section 
the corresponding identity for the excursion bridges is discussed. Here we 
start with by relating the distribution of (—go, d ) to the Levy measure of the 
inverse local time preparing in this way the connection of the Ito excursion 
measure and the Palm measure. In Section l3~2l we make some observations 
concerning the spectral representations of d and do — go- The result concern- 
ing excursion bridges is proved in Section EP1 bv utilizing the connections of 
the Ito excursion measure and the Palm measure. In Section 0] we discuss 
the null recurrent case and, finally, in Section [HI we give some examples. 



d t -t 
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2 Stationary excursions straddling t 



2 . 1 Preliminaries 

Let X° = {X t ° : t > 0} be a one-dimensional recurrent conservative diffusion 
living in the interval /. We assume that I = [0,6], 6 < +oo, or I = [0,6), 
6 < +oo, where is a reflecting boundary and 6 is either reflecting or natural 
or entrance-not-exit. For background on one-dimensional diffusions we refer 

to [131 and IB] 

The diffusion X° is characterized in I by its scale function S(x) and speed 
measure m(dx). Suppose, moreover, that X° is positively recurrent, that is 
M := m{I} < oo. Let p(t;x,y) be the symmetric transition density of X° 
with respect to m(dx) and let 

d d 
dm dS 

be the infinitesimal generator of X. Assume that m(dx) = m(x)dx, S(dx) = 
S'(x)dx such that m(x) and S'(x) are continuous and positive. The Green 
function is defined as 

POO 

G a (x,y) := / e- at p(t;x,y)dt. 
Jo 

Let P x and E x denote the probability measure and the expectation, respec- 
tively, associated with X° started at x and H y (X°) denote the first hitting 
time of y for X°. Recall from [13J p. 129 that 



E ( -aH v (X°)\ _ G »(x,y) 

x[ } ~G a (y, y y [Ll) 

Let now {X t (1) : t > 0} and {xf ] : t > 0} be two copies of X° such that 
Xq 1 '* = Xq 2 '* with the common law fh(dx) := m(dx)/M but let X^ 1 ) and X^ 
be otherwise independent. Define for t G R, 



X t (1) , t>0, 
xi 2) , t < 0. 



The process X = {X t : t G R} is called a stationary diffusion in stationary 
state living in I and having the generator Q. Notice that the law of X t is m 
for every t G R. 
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2.2 Joint distribution of go and do 

For the process X defined above let go and do be as in (jl.ip . In this section we 
derive the joint distribution of —go and do and show that the pair (—go, do) is 
an element of a special class /C of two-dimensional random variables studied 
in and [Ej. We start with the definition and some properties of this 
class /C. 

Definition 2.1. A two-dimensional random variable (^1,^2), where both £1 
and £2 are non-negative, belongs to K if 

fab) ±(U,V-U), 

where V is an arbitrary non-negative random variable and the conditional 
distribution of U given V is uniform on (0, V). 

Proposition 2.2. (Properties of the elements in K.) 

1) Let (£1,^2) £ K,. Then £1 = £2, £1 + £2 = V, and the density of £1 (and 
£2) exists and is given by 



F^x)= I v- l F v {dv). 

J fx, 00) 



[x,oo) 

If the density F v (v) =: p(v) exists then 

P(f 1 G dx, £ 2 G dy) = P ^ X + ^ dx dy. (2.2) 

x + y ' 

2) Let £1 and ^2 be non-negative random variables. Then (£1,^2) £ if and 
only if for all a ^ (3, 



E ( e -«&-«») = _J_ r E ( e ^+&)) d 7 . 



(2.3) 



Proof: For the proof, see and [Ej. □ 

Proposition 2.3. TTie pair (— go^o) belongs to the class JC with the joint 
Laplace transform given by 

E fp-^o+^so^ f i i I (0 A) 

[ > M(a-(3)\G a (0,0) Gp(0,0)J- K ' ) 
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Moreover, the joint density of {—go, d ) exists and is given by 
P(d G dt, —g G ds) 

p(t + s;y 1 ,y 2 )) dtds, (2.5) 



2/1,3/2=0+ 



M \dS{ Vl ) dS(y 2 

where p(t; x, y) is the transition density of the process X obtained from X° 
by killing X° when it hits zero. 

Proof: Using conditional independence of X^> and X^ 2 \ formula (|2.1|) . and 
the Chapman-Kolmogorov equation, we write 

E ( e - ad0+ ^) = ^m(dx)E :c (e-^( x(1) ))E a; (e- /3Ho ( x(2) )) 
m(dx) G a {x, 0)G/3{x, 0) /{G a {0, 0)G P (0, 0)) 

POO 

rh(dx) / dte~ at p(t;x } 0) 

POO 

x / ds e-P s p(s; x,0)/(G Q (0, 0^(0,0)) 



jo 

OO POO 

dt / dse~ at ~ ps 

JO 

x Jm(dx) p{t- x, 0)p(s; x, 0)/(G Q (0, 0)G (3 (0, 0)) 

OO POO 

dt / dse- at - (3s p(t + s;0,0) /(G a (0, 0)^(0, 0)) 
o Jo 

1 1 



1 

M 



M(a-(3) VG a (0,0) G p (0,0), 

To show that (—go, do) G /C we use Proposition 12.21 Letting a — ► /3 in (|2.4J) . 
it is seen that the derivative 75-Qy ) exists and 



E ( -ado+Pgo\ = 1 /" d ( 1 

1 ' M(a-/3) ^ c/7VG 7 (0,0) 



c?7 



x-P)Jb 



M(a 

To compute the joint density recall the formula (see [13], P- 154) 

P X (H (X°) G di)/<ft = = : ^(t ; x , 0) (2.6) 

dS (y) y =o+ 
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and, hence, 

P(d G dt,—g G ds)/dtds 

fh(dx)p^ (t; x, 0)p Jr (s; x, 0) 
d d 
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M \dS( yi ) dS(y 2 ) 



m(dx) p(t; x, yi)p(s; x, y 2 ) 



^/l,3/2=o^ 



d 



d 



M \dS{ yi ) dS{y 2 



-p(t + s;y 1 ,y 2 ) 



2/1,2/2=0+ 



where it is used that p(t; x,y) is continuous in t, x, y and that p + (t; x,y) is 
continuous in t, x and right continuous in y (see [T3| p. 149). □ 

Remark 2.4. To find the density ()2.5|) we used formula ()2.6|) and the 
Chapman-Kolmogorov equation. There is also an alternative approach. For 
this, let G a (x,y) be the Green function for the killed process X. Then (cf. 
[HI p. 185) 



d 



d 



and we have 
1 



dS(x) dS(y) 



1 



G a (x,y) 



x,y=0+ 



G a (0,0) 



(2.7) 



1 



a-P \G a (0,0) Gp(0,0) 
1 d d 



a — P dS(x) dS(y) 
1 



a- P Jo 



^0 



dt (e 



Gp(x,y) - G a {x,y) 
d d 



d 



-at\ 
d 



dS(x) dS(y) 



x,y=0+ 

P{t;x,y)\ x 



dS(x) dS(y) 



p(t + s;x,y)\ 



x,y=0+ 



y=o+ 
dt ds. 



2.3 Main identity in law 

The first main result of the paper is presented in the following theorem. 

Theorem 2.5. Let X be a stationary diffusion in stationary state. Then for 
allt G K 3 the two-dimensional random variables (t — gt, d t — t) and (If, It ) 
are identically distributed and, from Proposition\ 



E(e 



) 



1 



1 



1 



M(a-P)\G a (0,0) Gp(p,0)J' 



(2.8) 
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P(4 +) G dr,7 t (_) G ds) 
ltd d 



M \dS( yi )dS(y 2 ) 



p(s + r; y u y 2/ 



drds. (2.9) 

2/1,2/2=0+ 



Proof: We prove the theorem only for the case when the state space is / = 
[0, +oo) and leave other cases to the reader. 

By the stationarity, it suffices to consider the case t — 0. To simplify the 
notation, let 1^ := 1^ and J(~) := Iq^ . For y G (0, oo) introduce 

/ , rH {X) nH (X) \ 

:= I exp ^ - a J l{o<x s < y }ds - (3 J l {Xs>y} ds) j . 

From the Feynman-Kac formula it follows that u(x) is the unique bounded 
smooth solution of the equation 



Qu{x) 



au(x), < x < y, 
f3 u(x), x > y 



with the boundary condition u(0) = 1. The function u(x), therefore, has the 
following form 



u(x) := 



A^ a (x) + <p a (x)/<p a (0), x < y, 
Btp p {x), x>y, 



for some constants A and B, where ip a and <p a are the increasing and decreas- 
ing fundamental solution, respectively, of the generalized differential equation 

Qu = au (2.10) 
such that (p a is bounded and the "killing" condition holds 

fe(0)=0. (2.11) 

Both u and its derivative u + should be continuous at y, i.e., 

A$ a (y) + (p a {y)/(pa{ti) = B$p(y) 

(2-12) 

A^(y) + ^(y)/^a(0) = B^(y). 
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From (l2~T2l . 

<Pa(o)($Z(v)<Pi3(v) -$a(y)<pp(y)) 

few 

$Z(y)<Pi3(y) - i>a(y)£p(y) ' 

where it is used that the Wronskian (a constant) 

i := - i>a(x)fit(x) (2.13) 

equals -i/'J (0)¥> Q (0). Consequently, at point y, 

u (V ) = „ . 

i>£{yWp(y) -^a{y)fip{y) 

The joint Laplace transform of and 1^ is given then as 

E(e- oJ( - , -^ (+) ) = / m(dy)(«(y)) 2 
Jo 



oo 



= (ft(o))'/ a W .(2.14) 

Let 

r (y) ■= 4>a(yWp(y) - $ a {y)<Pfi(y) 

denote the denominator in the expression for w(y) and notice that 

= (a-flMy^y), (2.15) 

where we used that m(dx) = m(x)dx and S(dx) = S'(x)dx and that ip a ( x ) 
and (pp{x) satisfy the equations 

Qi) a = oap a and Q(p p = Pfip, 
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respectively. The equality (|2.15j) is the key to obtain 

{<P(i(y)y 



i>oo 

E(e-^ ( - , -^ (+) ) = (V>(0)) 2 / m(dy)- 

Jo 



( r iy)) 



2 



M(a-/})J Mv){r(y)Y 

M(a-I3),^\l ' yl My)(r(y)) 2 

Let F(e) denote the integral above. Using that (pp{x) is bounded and ip a (x) 
+00 as x — > +00, we have 

= r™{dyf M d 



$ a (y) dm V Kf) 

m (dy)-^- (-^) rs{dx) d i^ {x) 



dm \r(y)J J y dS \$Jx) 

S(dx)4-(S0] r m (dy) d ( 1 



S(dx) 



dm \r(y) 
1 1 



Since 



we obtain 



r(x) = -fy a (x)) 



dS \ip a (x) J \r(x) r(e) 
2 d I (pp{x) 



dS \Mx) 



F(e) = -I 5 dX \+ M£) 



(4> a (x)) ip a {e)r{e) 
By definition (|2.13^ of the Wronskian, 

1 d I (p a {x) 



2 • 
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y«(g)^g(e) -^(g)^(g) 
w a ($+(e){p p (e) - ij a (e)cp+{e)) 

Consequently, using boundary condition ()2.11j) . 
E(e-^- } -^ +) ) = ( 5 (0))2 x limF( £ ) 

fe(Q)) 2 (y a (o)y?(o) -^(Q)^(o)) 

M(a-/3) ^+(0)^(0) 

i f£?(Q) &(Q)' 

M{a-P) \<p {O) (p a {0) 
Finally, the Green function of X can be represented as 

G a (x, y) = w~ x ^) a {x)ip a {y), x < y 

(see, e.g., [H| p. 19), where i/; a , ip a are the fundamental solutions of ()2.1()j) 
such that ip a is bounded and 1])+ (0) = (condition for the reflection). Since 
fa — fa, we have 

1 <p+(0) 0+(O) 



(2.16) 



G Q (0,0) p a (0) £ a (0) 

Thus, the right-hand sides of ()2.16|) and ()2.4|) are equal and the proof is 
complete. □ 
Although we do not have any probabilistic explanation for the equality 
in law in Theorem 12.51 (but for RBM with drift, see the next section), the 
following lemma (cf. [22], Proposition 3.11 p. 330) explains why the variables 
It \ t — g t , and d t —t all have the same expectation. 

Lemma 2.6. For a given s > the events {X s > X } and {do > s} are 
independent. 
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Proof: By the symmetry of the transition densities, 



P(X S >X ;s<d ) = ^ ^Ip^xi 1 ) > x ; s < H (X^)) 



1 



b 



— I m(dx) I m(dy)p(s;x,y) 



J x 

b m(dy) f v ^ \ 
m(dx) p(s;y,x) 



o M Jo 

= P(X s <X ;s<d ). (2.17) 
Noting that, by the reversibility and the stationarity, 

P(X S > X ) = P(X. S > X ) = P(X > X s ) = 1/2 
completes the proof. □ 
Corollary 2.7. The expectations of and do are equal. 

Proof: By (12171) . 

E (/(+>) =2E( j\ {Xs>Xo} ds} = 2^ P(X S >X ; s < d ) ds 

POO 

= / P(d > s)ds = E(d ). 
Jo 

Similarly, we can show that E (jH) = E(d ). □ 

2.4 Probabilistic explanation of the identity when X is 
a RBM with drift 

In the case when X is a stationary reflecting Brownian motion with drift 

— fi < we prove the fact that = 1^ = d by considering the local time 
process of an excursion instead of applying the Feynman-Kac formula as in 
the proof of Theorem 12.51 

Let {L(t, y) : t > 0, y > 0} denote the local time of {X° : s > 0} up to t 
at the level y with respect to the Lebesgue measure, i.e., we have 

I g(X°)ds= [ g(x)L(t,x)dx (2.18) 

J0 J E 
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for any non-negative Borel function g. The following Ray-Knight theorem 
(see |Bj pp. 90-91) describes the behaviour of the local time process L up 
to H (X). Let Z^" 1 ' 2 ^ denote the squared radial part of an n-dimensional 
Ornstein-Uhlenbeck process with parameter ji and recall that the generator 
of Z in ' 2 ^ is 

d 2 , d 

2z _ + {n _ 2 „ 2) _ 

Theorem 2.8. Conditionally on X = x, 

{L(H (X),y) :0<y<x} = {Zf 2 ^ : < y < x}, 

{L(H (X),x + y):y>0}£ {Zf 2 ^ : y > 0}, 

where the process Z^ ' 2 ^ is started from the position of Z^ 2 ' 2 ^ at time x but 
otherwise Z^ 2,2 ^ and Z^ ' 2 ^ are independent. 

For the stationary excursion {X t : g < t < d } straddling zero introduce 
its total local time process by 

L«(y) := L<»{H {X<»),y) + L^(H (X^),y), y > 0, 

where L^(H Q (X^),y) and L^(H (X^),y) are the local times at y up to 
the first hitting times of zero for X^ and X^ 2 \ respectively. Since X^ 1 ' and 
are independent, given X = x, it follows that and are also 
independent, given X Q = x. 

Recall (see [32] and [2] p. 72) that if and are two independent 
non-negative time-homogeneous diffusions then Y^ + Y^ is a diffusion if 
and only if the generators of Y^ and Y^ are of the forms 

d 2 ., . d „ n 

ax—^ + {bx + Ci)—, 2 = 1,2, 
dx z dx 

respectively, where a > 0, q > 0, b G R. The generator of Y^ + Y^ is 
From Theorem 12.81 we obtain now the Ray-Knight theorem for . 
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Theorem 2.9. Conditionally on X = x, 

{L^(y):0<y<x}^{Z^:0<y<x}, 

{LV(x + y):y>0}±{Z^:y>0}, 

where the process Z^ 2 ^> is started from and the process Z(°' 2 ^ is started 
from the position of Z^ 4 ' 2 ^ at time x but otherwise Z^' 2 ^ and Z^ ' 2 ^ are 
independent. 

Lemma 2.10. Let H (L) := mf{y > : L^(X + y) = 0}. Then the random 
variables H (L) and L^ e \X ) are exponentially distributed with parameters 
2/x and [i, respectively. 

Proof: Define 

M (1) := sup{X t : < t < d } and M (2) := sup{X t : g < t < 0}. 
The distribution of M W (i = 1, 2) is given by (see jH] p. 14) 

P,(M'''<,) = P^cW< ff , W ) = ||^||. 

Let MW := M« - x. Then 

P X (M® <y) = P,(M« < x + y) = ^ (x+y) _ { , 

where we used that for RBM with drift 

1 



Hence, by integration, 



poo 

P(H (L)<y) = / 2/ie- 2 ^P x (M« < y)P t {M® <y)dx 
Jo 

2^ ^ , dx 



e 2ii{x+y) _ I 



1-e 



-2/12/ 
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Next we show that L^ e '(X ) ~ Exp(n), that is L^ e \X ) is exponentially 
distributed with parameter jj,. Recall that for Brownian motion with drift 
— /i killed at zero 

G (x,x) = ±-(e 2 ^-l). 
The Laplace transform of L(H , x) is (see |6J p. 32) 

E r e -aL{H {X),^ = ^ 1 = P_ (2 lg) 

V ' G (x,x)m(x)a + 1 (l-e~ 2 ^)a + /i 



From (l2~T9l) we find 



E ( e -«iW(Xo)) = [°° 2 -2^f \ 2 rfx 

v io V(l - e- 2 ^)a + 



/i + a 

□ 



Proposition 2.11. T/ie random variables E + \ E ' , and do have the same 
law. 

Proof: We show first that = /(">. By Proposition EH and (15131) . 



00 



-{X 3 >X }t 



E(e' a/(+) ) = E(exp(-a J 1 

= E(exp(-a^ ? ^)^)), 



where ~ #cp(jz) and C ~ Exp(2/i) is the life time of Z(°> 2 ^. Note 

that (cf. p. 90) 

{Z (0,2^) (C _ y) . < y < C} i {Z (4,2M) (y) . < „ < r}) 

where r ~ Exp(2/j,) is independent of Z^ 4 ' 2m \ Consequently, by Proposition 

E3H 



E (e"" /(+) ) = E ( exp ( - a jf Z^ 2 ^ 

= E(exp(-a^ ? 4°: 2 ; ) ) ^ 

= E^exp^-aj Zf^dy^j 

= E(e- a/(_) ). 



15 



To deduce that = do, we adopt the argument in [SI] used in the proof 
of the Biane Imhof identity. First recall that 

Z s ds, 

where the process Z = Z^ 2 ^ is the solution of the SDE 

dZ t = 2^/z~ t dW t - 2fiZ t dt, Z ~ Exp(fi). (2.20) 

Introduce next 

A t := / Z s ds 



o 



and let at be the right-continuous inverse of A t . Since the quadratic variation 

lo 



of the local martingale Y t = f* \fZ~ s dW s is given by 



(Y,Y) t = [ Z s ds = A t , 
Jo 



it follows from Levy's characterization theorem that Y at is a Brownian mo- 
tion, say B t , started at zero and stopped at A^. Hence, for t < Aoo, 



Z at -Z = 2 ^Z s dW s - 2/W Z s ds 
Jo Jo 
= 2Y at — 2fiA at 

= 2B t -2fit. 

Letting t — > A^ gives Z /2 = , where b[^ := —B t + fit. Thus, taking 
into account that 

< Z at = Z + 2B t - 2^t, < t < Aoo, 

yields 

A^ = M{t : = Z /2}. (2.21) 
Since Z ~ Exp(fi) and is independent of B^\ ()2.2H) is equivalent with 

/+ = 4» = H X (B^), (2.22) 

where A ~ Exp(2fi) is independent of B«. Noting that the right-hand side 
of (J2.22j) is identical in law to H (X) gives — do, as claimed. □ 
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3 Relation to Ito excursion theory 

3.1 Distribution of (—go, do) in terms of Levy measure 



Consider the distribution of do (and —go)- From (J2I 

E(e -* )=E( ^ H _1_ p.!) 

Next let £(t, 0) be the local time at zero up to time t (with respect to the 
speed measure) of {X s : s > 0}, X = 0. Let A = {A t : t > 0} be the 
right-continuous inverse of I. The process A is a subordinator and 

E (exp(-aA t )) = exp(-i *(a)), 

where the Laplace exponent $ is given by 

#(«)=/ (1- e- at )n + (dt) = a e~ at n + {t, oo) dt 
Jo Jo 

with the associated Levy measure n + (dt). 

Proposition 3.1. TTze Laplace exponent \l/(a) o/A «s giwen by 



y( a ) = - = - — E (e~ aH °) 



(3.2) 



x=0+ 

Proof: See p. 214. □ 
Proposition 3.2. Zet n + (dt) be the Levy measure of A. Then 

P(-9o e dt) = P(d e dt) = ra+( ^ r oo) dt; (3.3) 
v 

~P(V G (if) = — n + (dv), where V := d - g Q ; (3.4) 

i f 00 

P(d > v, -g > w) = — I n + (t,oo)dt; (3.5) 

M J v+W 

P(d G dt, —g G ds)/dtds = ~J^~^- n+ ft + s > °°)- (3-6) 
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Proof: Formula ([3.3)1 is a consequence of ([3.1)1 and Proposition 13.11 and is 
given in [23J and [2U ESI (the global formula). Formulae ([3.4)1 . ([3.5)1 . and 
()3.6|) follow immediately from (|3.3|l and Proposition 12.21 □ 



Remark 3.3. It is interesting to compare the right-hand sides of ()3.6|) and 
(pHj) . For this, notice first from El) that 



oo) 



d 



dS(x) 



P*(#o > *) 



a;=0+ 



(3.7) 



Thus, we have 
d 



— — n + (t + s, oo) 

Lib 



d d 
dt dS(x 

d d 
dS(x) dt 

d d 
dS(x) dS(y) 



P x {H <t + s) 
P x (H <t + s) 



2=0+ 



x=0+ 



p(t + s; x, y) 



x,j/=0+ 



3.2 Spectral representations for do and V = do — go 

In this section we show that the distribution of do (and —go) is a mixture of 
exponential distributions and the distribution of V = d — go is a mixture of 
gamma distributions. The mixing measures are the same and closely related 
to the so called principal spectral measure of the process X, as defined in 
Krein's theory of strings, see [TJJ [TH [THJ • Our starting point is the result due 
to Knight (see [T5] ) which states that the Levy measure n + (dt) is absolutely 
continuous with respect to the Lebesgue measure and there exists a unique 
measure A such that 



POO 

u(t) :=n + (dt)/dt= e- zt A(dz). 
Jo 

Moreover, A has the properties 



and 



A(dz) 
o z(z + l) 



A{dz) 



< oo 



oo. 



(3.8) 



(3.9) 



(3.10) 
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We remark (cf. jE]) that (|3.9|) is equivalent with the defining property of 
the Levy measure of a subordinator, i.e., 

(1 M)n + (dt) < oo. 

The property in ()3.10|) is a consequence of the fact that the speed measure is 
strictly positive everywhere and, in particular, in a right neighbourhood of 
0, see [Hj p. 82. For another proof of ()3.10|) see jH|, where A is interpreted 
as the principal spectral measure of a killed string. 

Proposition 3.4. Let A be the measure introduced above and associated with 
the process X.Then the measure 

A(dz) = A{dz)/{Mz 2 ) 

is a probability measure. Moreover, 

POO 

Y{d Q edt)/dt= I ze~ zt A(dz), (3.11) 
Jo 

and 

poo 

P(Vedv)/dv= / z 2 ve~ zv A(dz). (3.12) 



Proof: Recall that in the recurrent case (see [2H] P- 220 and §Bj p. 20), 

lim aG a (x, x) = — r— -, for all x G I. (3.13) 
a\o v ' m{I} v ' 

By (JUm, IpTHjl . and Fubini's theorem, 

1 

M := m{I} = lim 



a\o aG a (0, 0) 
n + (t, oo) dt 

POO 

dt / v(s) ds 



dt / A(dz 
Jo 

A{dz) 



e ~zt 



z 2 



and, therefore, A is a probability measure. Formulae (j3.11|) and ()3.12|l follow 
now from Proposition I3.2I and spectral representation (j3.8j) . □ 
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Remark 3.5. From the proof of Proposition 13.41 a new test for positive 
recurrence emerges: a recurrent diffusion X is positively recurrent if and 
only if 

A(dz) 

— V" 1 < oo. 
3.3 Excursion bridges 

In this section we use the theory of the Palm measures to show that the 
result in Theorem 12 .51 has a counterpart in the framework of excursions from 
to of the stationary non-negative diffusion X. Let 

M := {t G R : X t = 0} 

be the zero set of X and 

C:={teM : 3e>0V0<s<£ X t+S > 0} 

be the set of the starting times of excursions from to on a path of X. 

We introduce next the space E consisting of continuous functions e : 
R + i — > R + such that e(0) = and for which there exists C = C( e ) > with 
the property e(t) > for < t < ( and e(t) = for t > (. The space (E, £), 
where £ is the cr-algebra generated by the cylinder sets in the usual way, is 
called the canonical excursion space for excursions from to (associated 
with X). For t e C define X^ = {X^ x ' t] : s > 0} where 

x (ex,t) . = i x t+s, fort + s<R 
[0, fort + s>R, 

and R := M{u > t : X u = 0}. Clearly, X^ G E and ((X^) = R. 

We follow now j2Hj and define the concept of equilibrium excursion mea- 
sure which is a particular case of the Palm measure construction. See (2H] 
for general results on Palm measures and further references. For background 
on Palm measures, especially in queueing, see, e.g., [T]. 

Definition 3.6. For B E £ let 

Q(B) := E (\{t : < t < l,t e C,X^ G B}\) , (3.14) 

where | ■ | denotes the number of elements in the set under the consideration. 
The measure Q is called the equilibrium excursion measure for the excursions 
of X from to 0. 
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The next proposition and its corollary are adopted from [2H] (Theorem 
p. 290) where a more general statement originating from [2~T| and [22] is 
given. Because of the key importance of this result for our application, we 
restate and formulate it especially for excursions (cf. Section III in [2H] where 
connections with the Maisonneuve formula are discussed). 

Proposition 3.7. The measure Q is a-finite and for all measurable f : 
Rx E ^ [0,oo), 

e(^/(*,X^>)) = ! ! dsQ(de)f(s,e). (3.15) 
tec J R ^ E 

Proof: For the proof of the first claim, see [23]. To prove formula ()3.15|) . 
observe from the definition of Q, that ()3.15|) is equivalent with 




(3.16) 



Let 6 = {6t : t G M} be the usual shift operator defined in the underlying 
probability space via 

X a oe t (u)=Xt+ a (u) Vt,sGR. 

Notice that the law of X is invariant under 9; this is, in our case, equivalent 
with the stationarity of X. Consequently, by stationarity, 




= /" 1 ^e(^/( s ,x^)) 

^° sec 

= e(^/( s ,x^))), 

sec 
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hence (|3.16|) is proved giving (|3.15j) . □ 
Corollary 3.8. For a jointly measurable /i:Rx£h[0, oo) 

E (h(-g , X&'ri)) = [ Q(de) ( f * h(s, e) ds) . (3.17) 

J E JO 

In particular, for a > 0, v > 0, and e £ E 

P(-g e da, X (ex ' 9o) G de) = da Q(de) l( a <c(e)), (3.18) 

p( X ( e ^°) g de) = Q(tfe)C(e), (3.19) 
P(-3 Grfa)/da = Q((>a), (3.20) 

anc? 

P(y g <if) = v Q(C(e) e du). (3.21) 
Proof: To prove formula (I3.17J) . we set in (13. 15 J) 

/(t,e) = /i(-*,e)l(o<_t<f( e )) 

with /i jointly measurable (cf. |23j p. 291 and [22j p. 332). Then, due to the 
particular form of /, 

E(^/(t,X^)) =E (h(-g ,X^)) 
tec 

= / Q{de)( / h(s,e)ds\ (3.22) 

J E JO 

Choosing h appropriately in (l3~T7ll yields (l3~T8l . (EHTfll . and (l3~20l) . Finally, 
take in (|3.22D /i(s, e) = ls(e), 5 = {C > v }, u > 0, to obtain 

P(y>v) = Q(C(e);C(e)>«) 

which is equivalent with (|3.21j) . □ 

Remark 3.9. 1) From (l3~T8l) and (ITTTil) it follows that the distribution of 
— go given the excursion is uniform on (0, V) and, in particular, (—go, do) G /C 
(cf. Proposition 12. 3j) . 

2) It is seen from Proposition 13.21 and Corollary 13.81 that, 

Q(C(e) € *) = 
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To proceed, let M denote the Ito excursion law for the excursions of X 
from to 0. Then it is well-known that 



is the density of the first hitting time H := inf{t : X t = 0} given that 
X = x. Now (I3~271 yields 

M(C > a) = c'Q(C > a) 

for some (normalizing) constant c'. In fact, as can be deduced from (|3.18|) . 
the measures M and Q are the same (up to a constant): 

Proposition 3.10. There exists a constant d such that for all B G S 



We refer to [23j for a discussion about (|3.23j) and the associated normal- 
izations. See also [S] for the Brownian motion case. 

To formulate the main result of this section, we need the concept of 
excursion bridge of X. For reflecting Brownian motion the excursion bridge 
is a 3-dimensional Bessel bridge from to (of some length /). In general, 
the excursion bridge of X from to of the length / can be described as 
the process X( ''' ) obtained from X with X = x > conditioned to hit 
for the first time at time I by letting x —>■ 0. Clearly, defining X^ ' 1 ' ^ = for 
t > I it holds X(°>*'°) G E. 

Theorem 3.11. Let X^ 1 ^ = {xf' Z,0) : t > 0} be the excursion bridge as 
defined above and let U be a uniformly on (0, /) distributed random variable 
independent of X^°' l '°\ Let 



M(C > a) = c 




where c is a normalizing constant and 



n x (0, a) := P(H G da) /da 



M(B) = c'Q(B). 



(3.23) 




and 




Then L"' + > and ) are identically and uniformly on (0, /) distributed ran 
dom variables. 
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Proof: We use formula (|3.17j) and let therein 



h(s,e)=F(((e))G(^ * l {et>es} dt) 1 



{s<ae)} 

with Borel-measurable functions F and G. Recall the notation V = d Q — g , 
and consider first the left-hand side of ()3.17|) with h as above: 

r-V 

E (h(-g ,X^)) = E\F{V) G( J l [x i rao)>x ^ o)} dt) l { - go<v} ) 

= v(F(V)G(£l {Xgo+t>XQ} dt)) 

= E (F{V) G(JW)) , (3.24) 

where, as before, 

I« := / l { x t >x }^- 
-'so 

For the right-hand side of (|3.17l) we have 
Q(de)[j h{s,e)dsj 
= Q{F(((e)) C(e) J G[J l {et>es} d? 



C(e). 

= -M(F(C(e))C(e)y o G(jf l {et >e s} dt) ^) , (3.25) 

where in the last step Proposition 13.101 is used. Next recall that the descrip- 
tion of the Ito measure via the lengths of the excursions (see, e.g., [2E| P- 
497, [6J p. 60, [2Z| P- 421) says that for al\Be£, 



POO 

M(B) = / M(C E dl)P {0 ' l ' 0) (B), 
Jo 



where p(°''' ) is the probability measure associated with the excursion bridge 
process X^ ' 1 ' ^ defined in the canonical excursion space E. Consequently, we 
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obtain 



, K{e) Me) . J 

M(F(C(e)) C(e) J G[J^ l {et>es] dt) — 

M(C(e) G E^°)( j( G[ J l {et>es} dt 



ds\ 
TJ 



= J™ M(C(e) G dl)F(l)lE^(G^j\ {et>eu} dtj), (3.26) 

where {7 is uniformly on (0, 1) distributed random variable independent of 
X(°'''°). By lpTT7jl . the right-hand sides of f3~24]l and f3~25|) are equal. Con- 
sequently, using (|3.21|) in Corollary 13.81 and (|3.26j) . we obtain 

poo /»oo 

/ P(V G dZ)F(Z) E (G(/ (+) )|y = l)= P(V G dZ)F(Z) 
giving 

E = I) = (G(I^ +) )) . 

Combining this with the result in Theorem 12.51 completes the proof. □ 



4 Remarks on the null recurrent case 

Remark 4.1. Suppose now that X is null recurrent. Then the speed measure 
m(dx) serves still as the stationary measure of X but because m{I} = oo 
it cannot be normalized to a probability measure. However, the result in 
Theorem 12.51 is also valid in this case: 

Jm(dx) E^e-*™-^) 



In particular, given that V := do — go = v, 1^ and /( ) are uniformly 
distributed on (0, v). In the case when X is a stationary reflecting Brownian 



m 



(dx) E x (e- ado+ ^ 90 ) 

(4.1) 



a-p\GJ0,0) Gfl(0,0) 
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motion living above zero (|4.1|) gives 



^m(dx)E x (e- al(+) -^~ } ) = J m(dx)-E x (e- ado+ ^°) 

= -L-fy/to-y/W 
a — p V 



V2 



It is also easily checked that finiteness of m is not needed in the proof 
of Theorem 13.111 and, hence, for all diffusion excursion bridges X( 0,/ '°) of 
length I the random variables I^ + > and I^ 1 ^ (defined as in Theorem 13 .lip 
are uniformly distributed on (0, /). 

Remark 4.2. For a Brownian excursion there is an alternative simple proof 
of the result in Theorem 13.111 based on Vervaat's path transformation (see 
[SSHHED- Indeed, let X( 0,/ '°) denote a standard Brownian excursion of length 
I, that is a 3-dimensional Bessel bridge of length /, and U be a random variable 
uniformly distributed on (0,/) and independent of X^ ' 1 ' ^. Then the process 
X br = { X br - o < t < /} defined as 



x (o,i,o) _ x (o,i,o) ^ t + U < l; 
X^-X^, t + U>l 



equals in distribution to a standard Brownian bridge of length /. Clearly, 
with the notation in Theorem I3.11[ 

J0.+> = t l {xr>0} ds =: 



and 



o 



/ l {xl r <0} ds=:I^-\ 



o 



The claim follows now from the well-known result due to Levy, saying that 
j{br,+) an( j j(br-) are um f orm iy distributed on (0,/) (see, e.g., |6J p. 163 and 

[33 p. 43). 
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5 Examples 



To illustrate the results in Section [2] we consider some examples of stationary 
diffusions in stationary state. 

Example 5.1. (Reflecting Brownian motion with drift) Let X be a 

stationary reflecting Brownian motion with drift —fi < living above zero. 
The speed measure is 

m(dx) = 2e~ 2 ^ dx, 

the scale function is 

and the Green function at (0, 0) is 

G a (0,0) 



\j2a + fi 2 — /j, 
Hence by Theorem EH and (O (cf. [211), 

= ^((V2^TH?-/i)-(v^+A?-A*)) 
1\i 



^2a + /i 2 + y/20 + fi 2 
The transition density for the killed process is 

-vj. \ 1 u(x+v)-^— ( - {x ~ v)2 _ (*+v) 2 
p(t;x,y) = — ^— e ^ x +y> a e at — e 2* 

Thus by ()2.5|) . the joint density of (— go, do) and (j( + ), /(")) is 



= V c ~£(t+s) 



x,y=0 



(5.1) 
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To check the obtained formulae we give here also the following expressions 
for the corresponding Laplace transforms when X started at x: 

E x I exp ^ - a J l{o<x s <x} ds - f3 J l{x s >*} ds 

y/2a + y?e» x 



^2a + y 2 cosh(x y/2a + y 2 ) + ^2f3 + y? sinh(xy/2a + y?) 
(see [6J, formula 2.2.6.1 p. 300), while 



Example 5.2. (Brownian motion reflected at and 1) Consider now a 
stationary Brownian motion X living in the interval / = [0, 1] and reflected 
at and at 1. The speed measure of X is m(dx) = 2dx, the scale function 
S(x) = x, and the Green function at (0,0) is (see |6J, p. 122) 



'2a 

Hence the Laplace transform of (— go, do) and is (a ^ (3) 

v/^72tanh(\/2^) - ^/JT/2 tanh( y/2~j3) 



a — {3 



and 



E(e -a(do-3oh = — C s i n h(2v / 2^) + 2^ . 

4 v / 2^cosh 2 ( v / 2^) V 



Example 5.3. (Squared radial Ornstein-Uhlenbeck process) Let X = 

{X t : it 6 1} be a stationary squared radial Ornstein-Uhlenbeck process with 
parameters v = n/2 — 1 and 7. The generator of {X t : t > 0} is 

d 2 , n v d 
Q = 2x -— + {n- 27 x)—. 

dx z ax 

The process X can be defined as 

X(t) = e^SQBESie 2 ^^), t e R, (5.2) 
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where SQBES is a squared Bessel process of dimension n = 2u + 2 (see |24j). 
Assume that — 1 < v < and is a reflecting boundary. Then the process is 
positively recurrent with the speed measure given by (see [6] p. 140) 

m(dx) = -x u e~' yx dx, 

that is (after normalization) the gamma-density with parameters 7 and u+1. 
The Green function at (0, 0) is given by (see |H1 p. 141) 



G o (0,0) = 7 



where T(x) and B(x,y) denote the gamma and beta functions, respectively. 
By Theorem E3 

-aI<.+)-/3l(.-)\ _ -rp / - a d +i3g \ 



M(a -(3) VG a (0,0) (^(0,0) 
27 / 1 1 



a - /3 VB(§, -1/) B(£,-i/) 



Letting /3 = and taking the inverse Laplace transform (see [H] p. 261) 
gives the density of do (and 7^) and the expression for n + (t,oo)/M: 

P(4 6 *)/* = ^ = n _ v * 1 + v) *« (1 - e-"T • (5-3) 

Formula (|5.3jl can be also obtained using ()5.2|) as in [24J. The joint density 
of d and g ( an d 1^) is given then by 

f( t s ) = -L^L n +(t + s 00) = e ^(s+t) d _ e -2»(,+t)y - 1 
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